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We study the effect of electron interaction on the spin-splitting and the p-factor in graphene in 
perpendicular magnetic field using the Hartree and Hubbard approximations within the Thomas- 
Fermi model. We found that the (/-factor is enhanced in comparison to its free electron value g = 2 
and oscillates as a function of the filling factor v in the range 2 < g* < 4 reaching maxima at 
even v and minima at odd v. We outline the role of charged impurities in the substrate, which are 
shown to suppress the oscillations of the (/"-factor. This effect becomes especially pronounced with 
the increase of the impurity concentration, when the effective ^-factor becomes independent of the 
filling factor reaching a value of g* « 2.3. A relation to the recent experiment is discussed. 

PACS numbers: 72.80.Vp, 71.70.Di 



I. INTRODUCTION 

Graphene being subjected to a perpendicular magnetic 
field exhibits the unusual quantization of the energy spec- 
trum, which is manifested in a non-equally spaced se- 
quence of the Landau levels^. In contrast to conventional 
two-dimensional electron gas (2DEG) systems, the en- 
ergy difference between the lowest Landau levels is large 
enough allowing observation of the quantum Hall plateau 
even at room temperatures^. Another interesting pecu- 
liarity of graphene is the existence of the O'th Landau 
level located precisely at the Dirac point and equally 
shared by electrons and holesi. If the magnetic field is 
high enough, in addition to the Landau level quantiza- 
tion, the level splitting due to the Zeeman effect takes 
place. This kind of splitting was clearly observed in 
the recent experiments even for states lying relatively 
far from the Dirac point, at the filling factors v = ±4£ 
The Zeeman splitting is by its nature a one-electron ef- 
fect, which tells that a particle possessing a spin degree of 
freedom acquires the additional energy in the magnetic 
field B, 

V° z = cr 9f i B B, (1) 

where a = ±i describes two opposite spin states ^,-1; [is 
is the Bohr magneton, g is the free electron Lande factor 
(^-factor); g ~ 2 for graphene. However, experimentally 
observed splitting of the Landau levels can not be solely 
attributed to the Zeeman effect, as this splitting can 
also be enhanced by electron-electron interaction^. The 
electron-electron interaction in graphene is especially im- 
portant at high magnetic fields near v = when a new 
insulating state emerges^. Even though the nature of this 
state is still under debate, it is commonly believed that 
it is related to the electron-electron interaction^. 

The enhancement of the spin-splitting due the 
electron-electron interaction can be described by intro- 
ducing a phcnomenological effective (/-factor, g* , which 
effectively incorporates the interaction effects within the 
one-electron description. Calculation of the effective g- 



factor was originally done for conventional 2DEG sys- 
tems based on Si MOS^ and GaAs/AlGaAs^ structures. 
It was show that the (/-factor can be enhanced by the 
electron-electron interaction up to one order of magni- 
tude in comparison to its bare valued and oscillates as 
a function of a carrier density^. Interaction induced 
spin-splitting was extensively studied in confined 2DEG 
structures such as quantum wires^— . It was also ar- 
gued that interaction-induced spontaneous spin-splitting 
can take place in 2DEG systems even in the absence of 
magnetic fielcU^— . 

The enhancement of the effective (/-factor was also 
observed in carbon based systems. In graphite the ef- 
fective (/-factor is reported to be g* w 2.5^ Recently, 
Kurganova et. al£Q performed measurements of the ef- 
fective (/-factor in graphene. It was found to be g* = 
2.7 ± 0.2, which is larger than its non- interacting value 
g = 2. This indicates that electron-electron interaction 
effects play an important role and should be taken into 
account for explanation of the enhanced spin-splitting. 
Motivated by this experiment we use the Thomas-Fermi 
approach to study the spin-splitting in realistic two- 
dimensional graphene sheets in perpendicular magnetic 
field situated on a dielectric surface and subjected to 
a smooth confining potential due to charged impurities. 
(Note that the enhancement of effective (/-factor in ideal 
graphene nanoribbons has been recently studied by Ih- 
natsenka et. a}2i). The paper is organized as follows. 
Sec. II presents the model, where we specify the system 
at hand and define the Hamiltonian. In Sec. Ill we dis- 
cussed the obtained results and provide an explanation 
for the observed behavior of <?*-factor. Sec. IV contains 
the conclusions. 



II. MODEL 

We consider a system depicted in Fig. [TJ consisting 
of a graphene sheet located on an insulating substrate of 
the width d with the dielectric constant e r . (We choose 
e r = 3.9 corresponding to Si02). A metallic back gate 
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to Eq. @, the carrier density n a (r) at the position r 
depends on the total potential only at that position. 
The total potential 




V a (v) = V H (r) + VS(t) +V% + V mp (r) 



(4) 



is a sum of the Hartree, Hubbard, Zecman and the exter- 
nal potential produced by the impurities. The Hartree 
potential is given b y 26 ' 27 , 
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FIG. 1. (Color online) Schematic illustration of the studied 
structure. A graphene sheet is located on an insulating sub- 
strate of the width d separating it from a metallic gate. The 
substrate is contaminated by charged impurities with q = ±1 
situated at the distanace h — 1 nm apart from the graphene 
layer. 



is used to tune the carrier density by varying the gate 
voltage V g . We assume the charged impurities with the 
concentration n, are randomly distributed in the sub- 
strate at the distance h = 1 nm apart from the graphene 
layer—. The whole system is subjected to the perpendic- 
ular magnetic field B. In order to find the ground-state 
carrier density, we use the Thomas-Fermi approximation 
with the local relatior 
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v /|r-r'| 2 + 4d 2 



(5) 

where n(r) = n-j-(r) + n^(r) is the local carrier density, 
and the second term describes a contribution from the 
mirror charges^. The second term in Eq. (01 is the 
standard Hubbard potential which is shown to describe 
carbon electron systems in a good agreement with the 
first-principles calculation o 27 ! 29 



V£(r) = Urf (v)S a , 



(6) 



,23-25 



ra«r(r) = < 



/ p a {E-V{v))f FD {E-E F )dE (electrons) 

V"(t) 
V(t) 

J p a {E-V a {v))f^ D {E-E F )dE (holes) 



where U is the effective Hubbard constant and S a = 
is the area of unit cell of graphene(a ps 0.142 nm is 
the carbon-carbon distance). In our work we use U = 9.3 
eV which has been recently calculated within the con- 
strained random phase approximation^. The third term 
Eq. @ is the Zeeman energy given by Eq. (|T|). The 
last term in Eq. Q corresponds to the potential due to 
charged impurities and is given by 




^/\r-r. l \ 2 + (2d-hy 



(2) 

between the spin-dependent carrier density n CT (r) of the 
graphene and the total potential energy V{r). Here 

f FD (E-E F ) = 1/ (cxp (^) + l) and f h FD {E,p) = 

1— f' FD (E — E F ) are the Fermi-Dirac distribution func- 
tions for electrons and holes respectively, Ep = eV g is the 
Fermi energy. The Landau density of state in graphene 
is given by^ 



(7) 



where the summation is performed over charged impuri- 
ties in the dielectric; is the coordinate in the graphene 
plane of the projection of the i-th impurity of the charge 
qi situated at the distance h from the plane. Equations 
© and (Ql are solved self-consistently until a convergence 
is achieved. 

We define the effective g-factor as follows, 
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g*p B B = (vHv) - vHr)) 



(8) 



(3) 



where ui c = s/2vf/Ib is the cyclotron frequency, Ib = 
^Jh/eB is the magnetic length, vp is the Fermi velocity in 
graphene; the factor g v = 2 takes into account the valley 
degeneracy for all levels except of the zeroth one. The 
zeroth Landau level belongs both to electrons and holes 
which we take into account by setting g v = 1. According 



which assumes that spin-splitting in the system is caused 
by the Zeeman term, Eq. ([1}, where the free electron 
value g is replaced by the effective g- factor, g* . (If the 
Hubbard interaction is absent, U — 0, then apparently 
g* = g). Substituting Eq.(|4]) into Eq.®, we arrived at 
the equation used to calculate g*, 



(J* =9 + 



U 

fJL B B 



(n ; (r) -n T (r)) , 



(9) 



where (...) denotes spatial averaging over the graphene 
lattice sites. 



III. RESULTS AND DISCUSSION 

Figure presents the central result of the paper. It 
shows the effective (/-factor as a function of the filling 
factor v = — (riR = —73-, note that our definition of nn 
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includes a factor of 2 accounting for the valley degener- 
acy) calculated for different concentrations of impurities. 
The dependence g* = g*(v) exhibits two main features. 
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FIG. 2. (Color online) The effective g-factor as a function 
of the filling factor v for different concentrations of charged 
impurities, m = 0%, 0.02%, 0.08%, 0.2%, at the constant per- 
pendicular magnetic field B = 35T. Inset: the dependence of 
g* on the Hubbard constant U for the fixed v = 2.5. All the 
calculations are done at the temperature T — 4 K. 

First, the effective gr-factor is enhanced (g* > g) and os- 
cillates in the range 2 < g* < 4 achieving its maximal 
values at even filling factors v = 0, 2, 4, . . ., while having 
the minima at odd filling factors v = 1,3,.... Second, 
the increase of the impurity concentration suppresses the 
enhancement as well as the oscillatory behavior of g* , 
such that for high m the effective g-factor becomes only 
weakly dependent on v reaching the value g* « 2.3. 

In order to understand the observed behavior let us 
first consider in details the case of a low concentration of 
impurities shown on Fig. [31 where <?*-factor dependence 
for Hi = 0.02% is complemented by the spin-density and 
the polarization dependencies. For small rii the effect of 
impurities is small and the filling factor can be directly 
related to the number of the occupied Landau levels of an 
ideal system (i.e. without impurities). For a fixed value 
of the magnetic field the increase of the filling factor cor- 
responds to the increase of charge density through subse- 
quent population of the Landau levels. As seen in Fig. [31 

the total spin polarization P = (P(r)) = ( felTO^j} ) 

exhibits the same qualitative behavior as the <?-factor (ex- 
cept for v = 0, which will be discussed below). At v = 1 
the g* -factor reaches its minimal value g* = g. In this 
case the Fermi energy is located in between the 0'th Lan- 
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FIG. 3. (Color online) The dependence of (a) the charge con- 
centration, (b) the polarization and (c) the effective (j-factor 
on the filling factor in almost ideal system (i.e. in a graphene 
sheet with the low concentration of charged impurities in the 
substrate, n* = 0.02%). 



dau level (LL0) and the l'st Landau level (LL1), i.e. the 
LL0 is fully occupied, while LL1 is completely empty 
(see inset in Fig. [3|. This gives rise to the equal spin- up 
and spin-down densities and hence to the zero spin po- 
larization. The increase of the filling factor in the range 
1 < v < 2 leads to gradual population of l'st spin-down 
(1) Landau level (LL1(4)) and, in turn, to the increase 
of nj,, while does not change. (Note that even though 
in our model the DOS is given by the delta functions, 
it is effectively smeared out by a non-zero temperature, 
which results in a smooth change of the charge densi- 
ties.) Since the difference n± — increases, according 
to Eq. ^ g* grows and reaches its maximum g* w 3.5 
at v = 2, when the Fermi energy lies in the middle of 
two spin-split levels corresponding to the same Landau 
level (LL1). The enhancement of the effective g- factor 
in comparison to its noninteracting value is apparently 
caused by the Hubbard term in Eq. (0]). The Hubbard 
interaction enhances the spin-splitting triggered by the 
Zeeman interaction giving rise to g* > g. 
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When the filling factor is further increased from v = 2 
to v = 3, i.e. the Fermi energy is shifted towards higher 
energies, the population of the spin-up (f) level belonging 
to LL1 gradually grows, while the density of the spin- 
down electrons (\) belonging to the same LL1 remains 
unchanged as the later level remains completely filled. 
Eventually, at v = 3 the spin densities become equal, 
n i ~ ^f; the system is not spin-polarized (P = 0) and 
the effective (/-factor again reaches its minimum g* = g. 
The same physics is responsible for similar oscillatory 
behavior of the effective (/-factor and the polarization for 
higher filling factors. 
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FIG. 4. (Color online) Schematic illustration of Landau levels 
population at v = for electrons and holes. Shaded regions 
(red and blue) correspond to states occupied in the LLO by 
spin-up and spin-down electrons. 

The dependencies of the effective g-factor and the po- 
larization are qualitatively different for v < 1. Namely, 
the polarization drops to zero at v = 0, while the effective 
g- factor reaches its maximum, see Fig. [3jb),(c). This is 
in contrast to all other even filling factors when both g* 
and P exhibit maxima. This can be understood as fol- 
lows. In contrast to other Landau levels, LLO is equally 
shared by electrons and holes at Ep = 0, which is a dis- 
tinct feature of graphene. As illustrated in Fig. [H when 
the magnetic field is high enough, i.e. the spin-split lev- 
els are well resolved, electrons predominantly populate 
the LL0(|) state, while LLO(f) is mostly occupied by 
holes. As a result, = — n^, and therefore the effective 
(/-factor reaches the maximum because of the Hubbard 
term ~ U (nj, — 7if ) = 2U n^. On the other hand, at v = 
the graphene is electrically neutral, n — + = 0, 
and spin-polarization is absent, P = 0, since |rt^| = \rif\. 
Note that the effect of electron-electron interaction on 
spin splitting in graphene nanoribbons at v « was dis- 
cussed in Ref . . 

The above analysis is strictly speaking applicable only 
for ideal graphene, when m = 0. In this case, the range 
of (/-factor oscillations can be easily estimated from Eq. 
(J9]). At odd filling factors, when nj. = n-j-, Eq. Q gives 
9min = .9 = 2. At even filling factors, when Ep lies 
between two spin-split levels of a given Landau level, for 
the chosen parameters U and B, the effective (/-factor 
5max — 2 + _^ Sa _ « 4, which is in accordance with our 
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FIG. 5. (Color online) Distribution of (a),(e) the self- 
consistent potential, (b),(f) spin-up and (c),(g) spin-down 
electron density, (d),(h) the spin polarization for a single im- 
purity at different filing factors v — 2 and v = 3 (left and 
right columns respectively). 



numerical calculations (Fig. |2J n< = 0). However, in 
the presence of impurities, this is not the case anymore, 
as the oscillations of the effective (/-factor get suppressed 
and g* never reaches (/^ ax and always stays larger than 
2mi„> see Fig. [U 

In order to explain the influence of impurities on the g- 
factor, let us now consider a system consisting of a single 
repulsive impurity only. Figure [S] (a) shows the cross 
section of the self-consistent potentials and for 
spin-up and spin-down electrons respectively. The LL0(^ 
, f) coincides with the self-consistent potential V^^, while 
the positions of the LL1(|, t) are given by V^'^ + haj c . We 
have chosen two representative values of the filling factor, 
namely, v = 2 and v — 3 corresponding to maximum and 
minimum values of the effective (/-factor. 

At v = 2, which in ideal graphene corresponds to the 
almost occupied spin-down and almost empty spin-up 
states of the LL1, g* reaches the maximal value. Figure 
[5] shows that the LLl(^) is pinned to the Fermi energy 
Ep. (For the effect of pinning of Ep within the Landau 
levels see e.g. Ref4i). The states lying in the interval 
\E - E F \ < 2nk B T are partially filed < f FD < 1 and 
therefore the electron density can be redistributed un- 
der an influence of an external potential. These states 
represents the compressible strips^, which in our case 
extends over the whole system (except of the impurity 
region). The presence of negative impurity leads to the 
distortion of the potential as depicted in Fig. [5ja). As 
a result, in the impurity region the LL1(4_) raises above 
Ep and this state becomes depopulated, Fig. E)Jc). (Note 
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that LLl(f) is practically depopulated even in an absence 
of the impurity, Fig. [5^b)). As a result, the spin den- 
sity difference, n^ — n^, decreases in the impurity region, 
which apparently leads to the decrease of P and g* in 
comparison to ideal graphene, see Fig. [5Jd). 

On the other hand, the influence of the impurity is op- 
posite for odd v. At v = 3 the system is predominantly in 
a unpolarizcd state, which is manifested by the minimum 
of g* . However, the distortion of the potential due to the 
impurity gives rise to the formation of a compressible 
strip around the impurity, where Ep intersects the LL1. 
This is clearly seen in Fig. [S] (e) where the compressible 
strip corresponds to regions where the potential is flat 
because of the pinning to Ep within the energy window 
\E - E F \ < 2nk B T (where < f FD < 1). Because of 
the partial filling of the compressible strip, the electron 
density there can be easily redistributed there. As a re- 
sult, the Hubbard interaction pushes up and depopulates 
the LLl(t) while the LL1(4_) remains populated, see Figs. 
Eff),(g). This leads to a local spin polarization around 
the impurity as illustrated in Fig.([5])(h). therefore, the 
overall polarization is no longer zero, (-P(r)) > and, 
hence, the effective g-factor does not drop to the mini- 
mum value, remaining g* > g^ ia - 

Summarizing, the influence of a single impurity 
is twofold: when the system is predominantly spin- 
polarized (even v), the impurity decreases the average 
polarization and the effective g-factor by locally pushing- 
up the Landau levels and depopulating them; in the op- 
posite case of a predominantly non-polarized system (odd 
is), the impurity leads to the local formation of the spin- 
polarized compressible strips, which instead increases the 
average polarization and the effective ^-factor. 

Having understood the effect of a single impurity on 
the average polarization and the effective g-factor it is 
straightforward to generalize the obtained results for an 
arbitrary concentration of impurities. The higher con- 
centration n.i, the larger influence of impurities on the 
average value of the spin-polarization and the effective 
g-factor. As a result, an increase of the impurity con- 
centration leads to the suppression of the amplitude of 
oscillations as shown in Fig. [2 

Note that for a sufficiently large impurity concentra- 
tion (in our case rii = 0.2%), the oscillations of g* get 
practically suppressed and g* becomes rather indepen- 
dent on the filling factor, see Fig. O This effect can be 
understood from a comparison of two distinct cases of 
low and high impurity concentration, rii = 0.02% and 
iii = 0.2%, see Fig. [6] When the impurity concentration 
is low {rii = 0.02%, two left columns in Fig. [5]), the self- 
consistent potentials produced by different impurities do 
not overlap and the system can be treated as an assembly 
of independent impurities. (The potential is flat every- 
where besides narrow regions close to the impurities, see 
(a)-panels for m = 0.02% in Fig. H>|). At v = 2 the pres- 
ence of impurities decreases locally polarization (dips on 
(c)-panel) , while at v = 3 the local polarization increases 
(peaks on (c)-pancl). 



However, when the impurity concentration is high 
(rii = 0.2%, two right columns in Fig. [5]), the poten- 
tials produced by different impurities start to overlap and 
the analysis in terms of a single impurity is no longer 
justified. A given value of the filling factor can not be 
associated with a certain number of the Landau levels, 
since the potential is strongly distorted in comparison to 
the ideal case ((a)-pancls for rii = 0.2% in Fig. and 
therefore electrons occupy different Landau levels ((b) 
and (c)-panels for rii — 0.2% in Fig. [B]). In fact, the de- 
viations in the potential and densities from those of the 
ideal case become so significant, so the difference between 
the cases of v = 2 and v = 3 is practically washed out 
(c.f. two right columns in Fig. [5]). As a result, the av- 
erage value of the polarization and the effective g-factor 
becomes practically independent of the filling factor. 

In the model used in our calculation the enhancement 
of the ^-factor is caused by the Hubbard term in the po- 
tential, Eqs. © and ©• Let us briefly discuss how the 
calculated value of g* depends on the Hubbard constant 
U. While we used value U ~ 3.5*22, the current litera- 
ture reports various estimations of U in the range 0.5t < 
U < 2t)22r— where t « 2.7 eV is the hopping integral in 
the standard p-orbital tight-binding Hamiltoniani. We 
calculated the dependencies g* = g*(v) for different val- 
ues of the parameter U and found that the results show 
the same qualitative behaviour and the calculated value 
of g* scales linearly with U. This is illustrated in the inset 
to Fig. [2] which shows a dependence of the effective g- 
factor on the Hubbard constant for a representative value 
of 1/ = 2.5. 

Let us now discuss the relation of our findings to the 
recent experiment. Measurements done by Kurganova 
et aZ.— exhibit the enhancement of the effective spin- 
splitting leading to the effective g-factor g* = 2.7 ± 0.2. 
Also, the enhanced effective g-factor was found to be 
practically independent on v. Our calculations show that 
for low impurity concentrations, g* exhibits a pronounced 
oscillatory behavior in the range 2 < g* < 4, and it 
becomes rather independent of v for larger m reaching 
a saturated value g* ~ 2.3. Our calculations therefore 
strongly suggest that impurities always present in realis- 
tic samples play an essential role in suppressing the os- 
cillatory behavior of g*. Note that in real systems the 
oscillations of g* can be smoothed by a number of addi- 
tional factors. The measurements of Kurganova et. aJ2£ 
were performed in a tilted magnetic fields and at large 
filling factors v > 6. In this case the distance between 
the adjacent Landau levels is comparable to the Zeeman 
splitting which results in stronger overlap of the succes- 
sive Landau levels and eventually leads to an additional 
smearing of g*. Therefore our calculations motivate for 
further studies of the effective g-factor close to v = 0, 
where the oscillatory behavior of g* is expected to be 
more pronounced. Our finding also indicate that the os- 
cillatory behavior of the effective g-factor is expected to 
be more pronounced in suspended samples where the in- 
fluence of charged impurities will be much less important. 
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FIG. 6. (Color online). The spin resolved potential, densities and polarization for different concentration charged of 
impurities(?ii = 0.2%, 0.02%) and for different filling factors (y = 2,3). The one-dimensional plots of (a) V^(y), (b) the 
spin-up and (c) spin-down charge densities, and (d) the spin polarization P(y) as a functions of y for x = 50 nm. Dashed lines 
correspond to the ideal system (without impurities), (e) The 2D plot of the spatially resolved spin polarization P(x,y) in a 
graphene sheet. The system parameters are N x = 800, N y = 46f , d = 10 nm, B = 50 T. 



Finally, it is noteworthy that spin-splitting of in 
graphene^ and graphene quantum dota^ was also ex- 
perimentally studied in a parallel magnetic filed. It was 
concluded that in this case the effective <?-factor does dif- 
fer from its free-electron value. This can be explained by 
the fact that in the parallel field the Landau levels do not 
form and therefore the interaction induced enhancement 
of the g*-factor is small. 



showed that for low impurity concentration g* oscillates 
as a function of the filling factor v in the range from 
Smin = 2 to 5^ ax ~ 4 reaching maxima at even filling 
factors and minima at odd ones. Finally, we outlined the 
influence of impurities on the spin-splitting and demon- 
strated that the increase of the impurity concentration 
leads to the suppression of the oscillation amplitude and 
to a saturation of the the effective ^-factor around a value 
of g* w 2.3. 



IV. CONCLUSIONS 



In this work we employed the Thomas-Fermi approxi- 
mation in order to study the effective g-factor in graphene 
in the presence of a perpendicular magnetic field taking 
into account the effect of charged impurities in the sub- 
strate. We found that electron-electron interaction leads 
to the enhancement of the spin splitting, which is char- 
acterized by the increase of the effective g-factor. We 
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